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Abstract. Recently much work in Gamma-Ray Burst 
(GRB) studies was devoted to revealing the nature of out- 
burst mechanism and to looking for GRB progenitors. Sev- 
eral types of progenitors were proposed for GRBs. Most 
promising objects are coUapsars, compact object binaries, 
Helium star mergers and recently discussed supernovae. In 
this paper we consider four proposed binary star progen- 
itors of GRBs: double neutron star (NS-NS), black hole 
neutron star (BH-NS), black hole white dwarf (BH-WD) 
mergers and Helium star mergers (He-BH). Helium star 
mergers are a possible outcome of common envelope evo- 



cosmological distances there are still two major difBcul- 
ties in understanding this phenomenon. First, we do not 
understand fully the physics of outburst. Although sev- 
eral models were proposed they all yet have to meet some 
severe constraints imposed by observations (i.e. releasing 
energies of 10^^-10^'* ergs in timescales as short as 10~^ s 
in case of some GRBs). Second we don't know what are 
the astrophysical objects leading to gamma-ray bursts, i.e. 
what are their progenitors? 

Recently a black hole accretion disk model of GRB 



lution i)f a cumpacL ubjecL euLtiilng Lliy envtilupy uf a giant 



with a heliiiui cuie. 



Usiiig pupulaLluu syuLhesis we calciilatt; uuuibei uf Llic 
binary piugeuiLuiH and show LliaL BII-WD and Heliuui 
star 



outburst has been favored (Meszaros, 1999; Brown et al. 
2000; Fryer et al., 1999). Progenitors leading to this 
model include coUapsars ( Woosley, 1993| ; MacFadyen and 



Woosley, 1999|) and binary stars: Helium star mergers 



mrigeis duniinaLe pupulaLiuu uf the piupuscd biuaiy 
progeniLuis. Gumpaiisuu ufdisLiibuLiun uf diffuicnL biuaiy 
mergers around galaxies they are born in, with localiza- 
tion of GRB afterglows in their host galaxies shows that 
only Helium mergers may be responsible for GRBs with 



(Fryer and Woosley, 1998), double neutron stars (Ruffert 



et al., 1997; Meszaros and Rees, 1997), black hole neutron 



star (Lee and Kluzniak, 1995; Kluzniak and Lee, 1998) 
and black hole white dwarf systems (Fryer et al., 1999). 
Also recently the supernovae gamma-ray bursts connec- 



observ yd afLeiglows while it excluders NS-NS, BII-NS and — 2000l' I Chevalier 



BH-WB 



LhaL all 



systems as (jllli piugemtoib. Assummg 
GRBs come from Helium star mergers and comparing 
numbers of Helium star mergers with observed BATSE 
GRB rate let us derive upper limit on GRB collimation to 
be 



tion received much attention (Paczyhski, 1999; Woosley, 
200C), however there is still no clear ev- 



idence tlm^_th£se_two_j5h^ are intrinsically corre- 

lated ( Graziani et al., 1999 ). 

The present paper is an extension of our previous stud- 



ies (iBelczynski and Buhk, 1999| ; [Bulik et al., 1999| ; |BeL 



czynski et al., 2000) and we aim here at a discussion of 



Key words: stars: binaries; evolution — gamma rays: 
bursts — stars: neutron 



1. Introduction 

The last decade brought a great breakthrough in gamma- 
ray burst studies. Observations of GRB afterglows in X- 



ray, optical and radio wavelength domains (Costa et al.. 



1997; Groot et al., 19971:) led to identification of GRB 



host galaxies (Groot et al., 1997a) and measurements of 



their redshifts. This has solved the long standing problem 
of their origin. While we learned that GRBs come from 

Correspondence to: K. Belczynski, kabel@camk.edu.pl 



viability of the proposed GRB binary progenitors. 

One way of telling which group of proposed binaries 
might be responsible for GRBs is to predict their num- 
bers and compare them to the observed number of GRBs. 
Population synthesis is a powerful tool for predicting num- 
bers of binary populations although it suffers from many 
uncertainties as some parameters of single and binary evo- 
lution are poorly known. Moreover, population synthesis 
works well in predicting the relative numbers of events, 
while calculation of absolute numbers requires additional 
assumptions. Using population synthesis method we cal- 
culate production rates of Helium star mergers, double 
neutron stars, black hole neutron star and black hole white 
dwarf systems. We use BATSE detection limit of observed 
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number of GRBs to compare with our predicted numbers 
of the binary progenitors. 

Another way of discerning among the binary progeni- 
tors is to compare their merger site distributions around 
host galaxies with locahzation of GRBs within host galax- 
ies. Again this may be accomplished by population syn- 
thesis method to simulate a given binary population and 
then placing it in galactic gravitational potentials one may 
trace each binary until its components merge due to grav- 
itational wave energy losses. We perform such calculations 
for different galactic potentials and then compare the re- 
sults to the observations of GRB afterglows and their po- 
sitions within host galaxies. 

In this work we extend our previous studies to include 
two more proposed binary progenitors: black hole white 
dwarf binaries and Helium star mergers. Moreover, we im- 
prove our population synthesis code and for consistency we 
also present the updated results for two previously stud- 
ied types of proposed progenitors: double neutron star and 
black hole neutron star systems. 

Population synthesis method has already been ap- 
plied to study compact object binaries in context of GRB 
progenitors. However most of authors have concentrated 
only on double neutron stars and black hole neutron star 



Table 1. Localization of GRB afterglows in relation to 
the host galaxies 



Bloom et al., 1999; Belczynski and Bulik, 1999; Bulik 



systems (Narayan et al., 1991; Phinney, 1991; Tutukov 
and Yingelson, 1993; [Portcgies Zwart and Sprecuw, 1996 ; 



et al., ] 999; Belczynski ct al., 2000). Only one group (Fryer 
et al., 199£) presented calculations including all types of 
proposed binary progenitors and also coUapsars. Our cal- 
culations besides the conclusions we obtain, may serve as 
a direct chec k and comparison for the results presented by 
Fryer et al., (|l999|) . 

In section 2 we present observational data on the dis- 
tribution of GRBs around their host galaxies, in section 
3 we describe the population synthesis code used in this 
paper and formation scenario for Helium star mergers, in 
section 4 we discuss the results, and we present our con- 
clusions in section 5. 



2. Observations Distribution of GRBs 

The discovery of gamma-ray burst afterglows by the 
Beppo SAX satellite lead to identification of GRB host 
galaxies, and to localization of the GRB events in rela- 
tion to these galaxies. A list of GRBs with afterglows and 
their projected distances from the centers of host galaxies 
is shown in Table |^. 

From Table g we see that GRBs take place not far 
from the centers of their host galaxies and in the case 
of GRB970508 the offset is only 0.01" off the center of 
the host galaxy. For a review of recent observations see 



GRB 


redshift z 


Offset AG 


R± [kpc] 


970228 


0.695 


0.30" 


3.8 


970508 


0.835 


0.01" 


0.15 


971214 


3.42 


0.06" 


2.1 


980703 


0.966 


0.21" 


3.5 


980613 


1.096 


0.8" 


15 


990123 


1.61 


0.7" (?) 


16(?) < 5 


990510 


1.60 


< 0.08" 


< 2 


990712 


0.434 


0.24" 


1.4 



scribes only long GRBs as only for these bursts afterglows 
were so far observed. 



3. The model 

We use the population synthesis code described in detail 
in Belczynski 



Bulik ( |1999D . The code was modified to 
allow for evolution of low mass and intermediate mass 
stars. 

3.1. Modifications and description of population synthesis 
code 

Initial conditions. All binary star initial parameters are 
drawn from the same distributions, but we have changed 
the limits of binary components masses. As of now, we 
use more accurate description of single stellar evolution 
(see below); now we let the primary (more massive com- 
ponent at beginning) to have mass in range 8.6-100 Mq. 
The lower limit has been set in order to make sure that the 
primary will explode in a supernova explosion (unless it is 
stripped of mass in early mass transfer/loss episode) and 
turn to a neutron star or a black hole. The range of sec- 
ondary mass component is restricted now to 1.0-100 Mq. 
The lower limit here ensures that the star has a chance 
to produce a remnant in the Hubble time. This way we 
study population of binaries which have the most chances 
to produce the proposed binary star GRB progenitors. 
Single stellar evolution. To describe evolution of single 
star we use analytical formulae of Eggleton et al., (1989) 
and Tout et al., (1996) which are a new addition to our 



code. We follow a given star through different stages of 
its evolution: main sequence, Hertzsprung gap, red gi- 
ant branch, core helium burning, and asymptotic giant 
branch. The star may become a naked Helium star due 
to wind mass loss and then we follow its evolution until 
it cools down to become a white dwarf or if it is massive 
enough to explode as a supernova and become a neutron 
star or a black hole. 

We introduce one important change to the formulae of 



(Bloom et al., 2000). Moreover, the host galaxies are typ- 
ically small, irregular , with intense star formation pro- 
cesses ( Fruchtcr, 200C ). The data presented in Table |l de- 



Tout et al., (1996). In our calculations we assume that the 
mass of a compact object formed in a supernova explosion 
is half of the helium core mass of exploding star. This 
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results in the mass range of compact objects stretching up 
to ~ 20Mq , which can be compared to the maximum mass 
of compact objects of « 3Mq in the original prescription 
of Tout et al., (1996). We assume that compact objects 
below 3 M0 are neutron stars and over this limit black 
holes (for discussion see Belczyiiski et al., ( 2000| )). 



Depending on its initial mass any star will turn to be- 
come stellar remnant, i.e. a white dwarf, a neutron star or 
a black hole. During late stages of evolution (Hertzsprung 
gap through asymptotic giant branch) stars are allowed to 
lose mass via stellar wind. We follow Tout et al., (1996) 
and include standard Rcimcrs wind mass loss rate for- 



mula (Kudritzki and Rcimcrs, 1978) and for very massive 
and large stars wc include luminous blue variable strong 



winds (Humphreys and Davidson, 1994). It is important 
to note that Helium stars are not allowed to loose mass in 
our present code, and that leads to overestimations of the 
most massive compact object masses. Inclusions of Wolf- 
Rayet type winds for the naked Helium stars would de- 
crease our maximum compact mass down to « 12-15 M©. 
Such a change, however, would not change qualitatively 
conclusions of the present study. 

Binary evolution. Binary evolution may change entirely 
the evolution of any of its components. For components 
of close binary systems, which will interact during the 
course of evolution, we calculate mass loss/gain during 
mass transfer/loss events. We include common envelope 
evolution, quasi-dynamic mass transfer and hyper accre- 
tion onto compact objects during common envelope phases 
as in our previous studies ( Belczynski and Bulik, 1999| ). 
Components losing their Hydrogen-rich envelopes during 
giant stages become naked Helium stars, while those gain- 
ing mass (most often during they main sequence life) are 
rejuvenated. Rejuvenation and naked Heliu m sta r evolu- 
tion is treated as described by Tout et al., ( 1996| ). 

Any mass loss from the system (either through wind 
mass loss or during mass transfer/loss events) changes the 
binary orbit. Following evolution of binary and its compo- 



nents we also include tidal circularization ( Portegics Zwart 
and Vt rbunt, 1996| ) and magnetic breaking (Tout et al 
1996()~| 

During supernova explosion we follow precisely the or- 
bit evolution and we check if a system in which a super- 
nova explosion takes place survives the event. We calcu- 
late the mass of the newly formed compact object, the 
remaining mass is expelled from the system carrying off 
momentum, and a natal kick is added to the newly formed 
compact object (either a neutron star or a black hole) and 
then the new orbit is calculated. Systems which survive 
explosions receive additional center of mass velocity as an 
effect of the natal kick and mass loss from the outbursting 
component. Once the evolution of binary components is 
terminated and a system has survived mass transfer/loss 
events and supernova explosions we study populations of 
proposed binary GRB progenitors, that is NS-NS, BH- 
NS and BH-WD compact object binaries. These binaries 



evolve only due to gravitational wave energy loss which 
will cause the orbital separation decrease and finally lead 
to a merging event of two compact components and pos- 
sibly to a gamma-ray burst. For supernova explosions and 
gravitational wave energy loss we use same treatment and 
formulae as in Belczynski and Bulik, (1999). 



3.2. Helium star mergers 

In this subsection we will describe more specifically an evo- 
lutionary path which may lead to formation of a Helium 
star merger, a new feature of our code. 

The binary components may merge during mass trans- 
fer events, provided that the orbit is not too wide. As an 
example let us consider a binary, with an small enough 
orbital separation, so that the primary (the initially more 
massive component) during its expansion on the giant 
branch overfills its Roche lobe. Let us suppose that the 
mass ratio (secondary to primary) is small, so that the 
mass transfer will proceed on the dynamical timescale. 
The entire envelope of the giant primary will be lost, and 
it will become a naked Helium star. The secondary which 
is still on its main sequence will not have the time to accept 
any of the matter shed from the primary envelope, so it 
will survive the mass loss virtually unchanged. As a conse- 
quence of this mass loss event the binary orbit will shrink 
drastically. If there is enough orbital energy to expel the 
entire envelope of the primary then the system will survive 
(otherwise the helium core will merge with main sequence 
star, a case we are not interested in this study) and con- 
tinue its evolution. If the Helium star is massive enough 
it will undergo a supernova lb type explosion. Since the 
orbit is tight after the first mass transfer there is a good 
chance that the system will survive the supernova explo- 
sion. Thus the system now consists of a neutron star or 
a black hole and a main sequence star. As a consequence 
of mass loss from the system and randomly added kick 
to newly formed compact object the orbital separation in- 
creases and the orbit becomes eccentric. As the time goes 
on, the secondary, which is probably a more massive com- 
ponent now, will start its evolution up the giant branch. 
Once the secondary radius approaches the Roche lobe, 
tidal interaction circularizes the orbit and decreases the 
orbital separation. When the secondary overfills its Roche 
lobe the system goes through a similar phase as during the 
first mass transfer. If, as we have mentioned above, the 
secondary is the more massive component then the mass 
transfer proceeds again on dynamical time scale. In this 
mass transfer phase when the envelope of the secondary 
engulfs the system in a common envelope the compact ob- 
ject begins to spiral in toward the helium core of the giant. 
There is a chance now that the compact object will accrete 
some material from the envelope, as pointed by Bethe and 



Brown, (1995 ), and if it is a neutron star it may collapse to 



form a black hole. Since the orbit is already tight when the 
system enters the second mass transfer phase the orbital 
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energy available is not sufficient to expel the envelope of 
the secondary. The spiral in continues until the black hole 
enters Helium core of the giant. It disrupts tidally the He- 
lium core, swallowing at the beginning a part of it, but the 
remaining helium material will form a hot, rapidly rotat- 
ing disk around the black hole ( Fryer and Woosley, 1998| ) . 
This configuration, a Helium star merger, may lead to a 
gamma-ray burst in the black hole accretion disk GRB 



model (Fryer et al., 1999) 



4. Results and discussion 

4-.1. Relative production rates 

In Figure |l] we show the relative numbers of four differ- 
ent GRB progenitor types that merge within the Hubble 
time (15 Gyrs) as a function of the width of the distribu- 
tion from which we draw kick velocity a compact object 
receives in a supernova explosion. 

Two things are clearly seen; first the number of WD- 
BH binaries and Helium mergers (He-BH) is about the 
same and is more then an order of magnitude greater then 
the number of NS-NS and BH-NS binaries. It means that 
if GRBs originate from binary mergers then they mainly 
come from WD-BH binaries and/or Helium mergers. 

Second, the relative number of a given progenitor type 
falls off approximately exponentially with the kick veloc- 
ity. This is quite clear, the larger kick compact object re- 
ceives the larger chance that the system will be disrupted 
in supernova explosion. And this is a reason for smaller 
production rates with increasing kick velocity. Assump- 
tions about the kick velocity distribution plays an impor- 
tant role on results of compact object binaries population 
synthesis. We draw the kick velocities from a three dimen- 
sional Gaussian distribution of a given width which we 
treat as a parameter in our studies. There is a line of ev- 
idence coming form pulsars galactic velocity observations 
that the distribution of kicks is bimodal and consists of a 
weighted sum of two distributions: about 80% are inter- 
mediate kicks of about 200 km s~^ and the other 20% are 



high veloc ity kicks of about 700 km s ^ ( Cordcs and Cher- 
noff, 1997). However if real distribution is in fact bimodal. 



then the high velocity component won't have significance 
for properties of the population of compact object bina- 
ries, as these will tend to be disrupted by high velocity 
kicks. High velocity component will imprint its presence 
in the kick distribution through high peculiar velocities of 
single pulsars and will decrease number of compact ob- 
ject binaries. The lesson from this is that if one studies 
properties of compact object binaries, and not their num- 
bers, then only the lower component of the kick velocity 
distribution is relevant. 

Relative production rates may be calibrated to obtain 
the real rat es in our Galaxy(eg. see eq. 14 in Belczynski 
and Bulik, (119991 )). For example, for the width of kick ve- 



0,1 



10- 



10- 



^10- 



10- 



lO-f ^ 




WD-BH 
He-BH 



NS-NS _ 



BH-NS 



200 



400 



600 



800 



1000 



Fig. 1. Relative fraction of different GRB progenitor 
types that merge (or manage to evolve) within the Hubble 
time. 



per Milky Way like galaxy per 10^ yrs for BH-NS sys- 
tems, 3 events for NS-NS binaries and 60 for WD-BH and 
Helium mergers. These numbers are obtained under as- 
sumption that binary fraction is 50%, that there are 0.02 
supernovae per year in the Galaxy, and that the star form- 
ing process has been constant throughout the history in 
the Milky Way. 

4-2. History of binary merging events 



The results of the population synthesis code can be com- 
bined with our knowledge of the star formation rate his- 
tory to yield the rate of various types of GRB progenitors 
as a function of redshift. Star formation history at high 
redshift is not well known, however it is generally agreed 
that the star formation rate rises steeply up to z w 1. 
At higher redshifts the analysis of the Hubble Deep Field 
( Madau et al., 1996| ) provided lower limits on the rate, yet 
these limits decrease with increasing redshift. On the other 
hand Rowan- Robinson ( 1999| ) argues that the star forma- 
tion does not fall down and remains roughly at the same 
level above z — 1. We consider two cases: a star formation 
function falling down steeply above z « 1 (the thin line in 
Figure H) , and a case of strong star formation continuing 
up to 2: = 10 (the thick line in Figure ^). 

For a given type i of the GRB progenitor we can cal- 
culate the number of events up to the redshift z: 



locity of Vki 



200 km s we obtain: 1 merging event 



Ni{> z) = 47r 



rf-^^dz. 

^ dz l + z 



(1) 
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Fig. 2. Star formation history rates used in this work. The 
thin hne is based on Madau et ah, (1996), while the thick 
line represents approximately the rate of Rowan- Robinson 
(1999). 

where is angular size distance — cHq + 
JIa)^/^. Ri{z) is the rate of a given type of event at the 
redshift of z: 



R,{z) 



t' 



Rsfr{z)*h*{t')p{t-t')dt\ 



where t 
0„)(14 



t(z) 

is the conformal time, 



(2) 



, dt = dz{l + z)-^((l + 
z)2 - z{z + 2)17a)^i/2^ p(^) is |.]^g probability dis- 
tribution of a merger of a given type as a function of time 
since formation of the system. We calculate the distribu- 
tion pit) for each type of a merger using the population 
synthesis method. 

In Figure ^ we show the cumulative rates of differ- 
ent merging events as a function of redshift. We have 
combined our relative numbers (shown on Figure H) for 
different progenitors with the star formation rate func- 
tion ( Madau et al., 1996| ) and ( Rowan-Robinson, 1999| ), 
and after taking into account the evolutionary time de- 
lay of a given merging event we integrated our relative 
production rates to find the merger rates as a function 
of redshift. In this example calculation we used two cos- 
mological models with and without the cosmological con- 
stant: = 0.3, JIa = 0.0, and VIm — 0.3, JIa = 0.7. In 
both cases the Hubble constant is Hq = 65 km s~^Mpc^^. 
We used the kicks drawn from the distribution which is a 
weighted sum of two Gaussians: 80 percent with the width 
of 200 km s""'^ and 20 percent with the width 800 km s^^. 

The curves in Figure ^ can be compared with the 
BATSE gamma-ray burst detection rate corrected for 
BATSE sky exposure, which is « 800 events per year. 
Comparison of the cumulative distributions for different 



progenitor types with the BATSE rate shows that if any 
of the progenitor types included in our calculations were 
to reproduce the BATSE rate, then we should not see 
GRBs from redshifts greater then z ~ 0.3! Of course this 
is not the case (e.g. see Table |l|), as GRBs with higher 
redshifts were observed, and the median observed GRB 
redshift is z ~ 1. However, we have not yet introduced the 
collimation factor into our results. The predicted cumula- 
tive rates presented in Figure ^ will decrease if we account 
for collimation and thus restricted visibility of gamma-ray 
bursts. To lower down our calculated rates to the BATSE 
rate, for average GRB redshift of about unity, we would 
need the collimation of about 4° for BH-WD and Helium 
mergers and about 12° for BH-NS and NS-NS mergers. 
And, as the population is dominated by BH-WD and He- 
lium mergers then overall requirement for collimation is to 
be « 4°. This is only an upper limit as we would expect 
such a collimation factor if all GRBs were originating only 
from binary progenitors. 

The thin lines in Figure ^ flatten out for high redshifts 
(z > 5). In other words we do not expect binary mergers 
at high redshifts. This is a combined effect of two factors. 
First, the star formation rate function (SFR) we have used 



(Madau et al., 1996; Totani, 1997) falls down steeply for 



high redshifts. This means that at high redshifts we do 
not expect many stars, and thus their mergers. However 
the SFR we have used is highly uncertain for high red- 
shifts and our result here may be quantitatively question- 
able. However, the thick lines in Figure^, corresponding to 
the predictions of the SFR of Rowan-Robinson (1999) do 
climb up with redshift. In this case the number of GRBs 
up to the redshift of 10 is nearly double that up to the 
redshift of 2. Thus, future detection (or non detection) of 
GRBs from such high redshifts will serve as a probe of the 
SFR at large redshifts. 

Non zero lifetimes of binary progenitors are the second 
thing that makes our curves to flatten out with redshift. 
Binary GRB progenitors need a specific time to evolve to a 
compact object binary or to a Helium merger {ievoi) and 
the compact object binaries need time to merge due to 
gravitational wave energy losses (tmerger)- This times are 
non negligible and are specific for each group of proposed 
binary GRBs progenitors. For our sample of binaries we 
found characteristic lifetimes which are the sum of the evo- 
lutionary times and merging times tufe [tevoi + ^merger)- 
They are, for NS-NS: ~ lO^-lO^^ yrs, for BH-NS: - 10^- 
10^" yrs, for He-BH: ~ lO^-lO^ yrs, for BH-WD: - 10^- 
10^^ yrs. We see that these times are non negligible, and 
even if star formation process has begun at some point, 
for a given z, we need to wait at least 10^-10^ yrs to start 
producing GRBs of binary origin. 

An additional point to emphasize here is that the evo- 
lutionary times {tevoi), besides merging times, may play 
important role for some types of binary GRB progenitors. 
They are less important for NS-NS and BH-NS systems 
which are end products of high mass stars (very fast evo- 
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Fig. 3. Cumulative event rates of different types of GRB progenitors. The left panel corresponds to the cosmological 
model described by = 0.3 and = 0.7, while the right panel to ilm = 0.3 and = 0. The solid line describes 
the NS-NS rate, the dash dotted line describes the BH-NS rate, the short dashed line corresponds to the BH-WD 
rate and the long dashed line shows the He-BH rate. The thin lines corresponds to the assumed star formation rate of 
Madau et al., (1996), while the thick lines represents the results obtained using the rate of Rowan-Robinson (1999). 
Note that BH-WD, and He-BH events strongly dominate the population. 



lution) and for which merging times {tmerger) Eire compa- 
rable with their total lifetimes tufe. However evolutionary 
times are important for BH-WD systems, as it may take 
a long evolutionary time to form a white dwarf, usually 
comparable with the merging time of final compact black 
hole white dwarf binary. It is even more clearly seen in 
the case of Helium mergers for which evolutionary times 
equal total lifetimes, as these systems merge during com- 
mon envelope evolution of still unevolved (not a compact 
object) binary. 

The above example calculation shows that the WD-BH 
and He-BH events are far more numerous than the NS-NS 
or BH-NS events. The absolute numbers presented above 
depend on a number of assumptions leading from the bi- 
nary population synthesis to the observed rate. The sec- 
ond rather robust result is that if the star formation rate 
remains high even at high redshifts, the number of GRBs 
from such high z must be significant - see the difference 
between thin and thick lines in Figure zdist. 

4-. 3. Distribution of binary mergers around host galaxies 

In Figure ^ we present the distribution of center of mass 
velocities gained by systems in the supernova explosions 
versus binary lifetimes (the time binary takes to evolve 
from ZAMS to final merger of two components) . Data is 
presented on four panels for four types of the binary GRB 



progenitors. On each panel we plot three lines to show the 
distribution of mergers in respect to host galaxy of the size 
and mass comparable to Milky Way. Together these lines 
define the region in the parameter space with systems that 
can escape from a massive galaxy ( Bloom et al., 1999 ). 

The horizontal dashed line corresponds to the Hubble 
time (15 Gyrs); we are not interested in the system above 
this line as they do not have a chance to merge within 
Hubble time. 

The vertical solid line corresponds to v — 200 km s~^, 
which is approximately the escape velocity from an in- 
termediate mass galaxy. Systems to the right of this line 
(with velocities higher then the escape velocity) will merge 
outside their host galaxy. 

The inclined solid line corresponds to a constant value 
of w X tmerge = 30 kpc, which is about the radius of a 
high mass galaxy. Systems above this line have gained high 
enough velocity and have enough time to escape from their 
host galaxy to merge outside of it. 

We note that a significant fraction of NS-NS and BH- 
WD merging events takes place outside of the host galax- 
ies, and thus their distributions are inconsistent with the 
GRB observations. This conclusion is even stronger be- 
cause it was drawn for the case of massive large host 
galaxy. As noted in Sec. 2, GRBs' host galaxies are small, 
and in this case even more NS-NS and BH-WD mergers 
would take place outside their hosts then it is inferred 
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Fig. 5. Cumulative distribution BH-NS events around 
low mass host galaxies. The solid line corresponds to the 
case of no natal kicks Oy — Okms^^, the dotted line shows 
the case oiay — 100km s~^, and the dashed line represents 
the case of = 200km s~^ More than 50% of BH-NS 
merger events take place outside the 15 kpc radius from 
the host. 
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Fig. 6. Cumulative distributions of Helium star merg- 
ers (He-BH) events around low mass host galaxies. The 
solid line corresponds to the case of no natal kicks — 
0kms~^, the dotted line shows the case of cr^ = lOOkms"^, 
and the dashed line represents the case of cr^ = 200km s"'^. 
Almost 95% of He-BH events takes place within the 15 kpc 
radius from the host. 



form Figure |^. One may argue that binary mergers tak- 
ing place outside galaxies in a thin intergalactic medium 



This excludes BH-NS systems as potential GRB progen- 



itors given that GRB host galaxies are small (Fruchter 



would not produce afterglows, and thus would not con- 2000) 



tribute to observed distribution of GRBs in relation to 



their host galaxies. However according to Costa (|2000|) all 
GRBs observed by Beppo SAX are accompanied by after- 
glows. 

Thus, we are left with two binary progenitors: BH-NS 
systems and Helium star mergers. As seen from Figure || 
we predict that all of them would produce gamma-ray 
bursts inside their host galaxies provided that hosts are 
large and massive. But as noted before this is probably 
not the case, and we have calculated their distribution 
in case when their hosts were small, low mass galaxies. 
We approximate trajectories of our systems in a poten- 
tial of a small galaxy by propagation in empty space. We 
place them in one point and let them move with veloci- 
ties gained during supernova explosions. We follow their 
trajectories until they finish their life in a merging event. 
We present the expected cumulative distributions of the 
projected distance from the center of the host galaxy for 
BH-NS merging events in Figure |^ and for Helium star 
mergers in Figure ||. 

In Figure ^ we see that more than 50% of BH-NS sys- 
tems merge outside the 15 kpc projected radius. This is 
in clear disagreement with the GRB distribution obser- 
vations, which show that all so far observed GRBs take 
place within 15 kpc radius off their host galaxy center. 



Thus the only possibility left for the binary origin of 
GRBs are Helium star mergers. From Figure ^ we see that 
almost 95% of these merging events take place within the 
15 kpc radius from the host, just as expected if they were 
gamma-ray burst progenitors! This result is independent 
of the host galaxy size and mass, as we have shown that 
these mergers will merge within 15 Kpc radius even in the 
case of no pull from their host - in the case of propagation 
in empty space. 

Now, if we assume that all GRBs are coming from bi- 
nary stars we must remark that all long GRBs result from 
Helium star mergers and we estimate collimation factor 
of these burst to be ~ 4°). All the other merging events, 
coming form NS-NS, BH-NS and BH-WD systems may be 
responsible for short GRBs, for which afterglows and thus 
their distributions in respect to host galaxies have not yet 
been observed. 

5. Conclusions 

We have calculated the properties of the possible binary 
GRB progenitors: BH-NS, NS-NS, BH-WD, and He merg- 
ers. The GRB binary progenitor production rates fall off 
exponentially with width of natal kick velocity distribu- 
tion. We calculate the expected redshift distributions and 
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Fig. 4. The distribution of different type of compact object binary mergers in the space spanned by the binary center 
of mass velocity and the binary hfetime (from Z AMS to final merger) . The horizontal dashed line corresponds to the 
Hubble time (15 Gyrs). In the region for tmerge < 15 Gyrs we present two solid lines: the vertical corresponding to 

V = 200 km s~^ - approximately the escape velocity from a galaxy, and the line corresponding to a constant value of 

V X tmerge = 30 kpc. Together these lines define the region in the parameter space with systems that can escape from 
the host large galaxy. 



numbers for each type of the progenitor, and find that 
BH-WD and He-BH type events dominate over the NS- 
NS or BH-NS mergers. Moreover, in calculating the red- 
shift distribution of binary progenitors of GRBs one can 
not neglect the evolutionary times tevoi- In our example 
calculation we find that assuming that all GRBs result 
from binary mergers, then the population is dominated 



by BH-WD and He-BH star mergers, and the coUimation 
must be of order 3 x 10"'^ (~ 4°). The existence or non 
existence of GRBs at high redshifts depends on the star 
formation history. In particular the recent measurement 
of z = 4.50 in case of GRB000131 ( [Andersen ct al., 2000D 
suggests that the star formation rate has been high up to 
high z. 
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We find that only the He-BH (Hehum star merger) 
model of GRBs is consistent with Beppo SAX observa- 
tions of long bursts and their localization in host galaxies. 
As noted in several previous papers a larger fraction of 
NS-NS binaries merge outside even massive host galaxies. 
The distribution of BH-NS mergers is tighter around the 
massive galaxies, however we find that they do escape from 
the potentials of small galaxies. Our population synthesis 
code results indicate that a significant fraction of the BH- 
WD binaries should obtain large enough velocities so that 
they can escape from the potential well even of a massiv e 
galaxy. This differs from the results of Fryer et al., ( 1999| ), 
who find that the BH-WD binaries should lie within the 
host galaxies. Moreover our code results in a much larger 
f ractio n of BH-WD binaries than the fraction Fryer et al., 
( |1999| ) obtain. Finally we remark that NS-NS, BH-NS and 
BH-WD mergers can be responsible for short bursts which 
are not observed by Beppo SAX. This could be resolved 
by HETE-H, if short burst afterglows are discovered, and 
if their host galaxies are identified. 
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